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Abstract 

"Einstein-aether" theory is a generahy covariant theory of gravity containing a dynamical pre- 
ferred frame. This article continues an examination of effects on the motion of binary pulsar 
systems in this theory, by incorporating effects due to strong fields in the vicinity of neutron star 
pulsars. These effects are included through an effective approach, by treating the compact bodies 
as point particles with nonstandard, velocity dependent interactions parametrized by dimension- 
less "sensitivities" . Effective post-Newtonian equations of motion for the bodies and the radiation 
damping rate are determined. More work is needed to calculate values of the sensitivities for a 
given fluid source; therefore, precise constraints on the theory's coupling constants cannot yet be 
stated. It is shown, however, that strong field effects will be negligible given current observational 
uncertainties if the dimensionless couplings are less than roughly 0.01 and two conditions that 
match the PPN parameters to those of pure general relativity are imposed. In this case, weak 
field results suffice. There then exists a one-parameter family of Einstein-aether theories with 
"small-enough" couplings that passes all current observational tests. No conclusion can be reached 
for larger couplings until the sensitivities for a given source can be calculated. 
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I. INTRODUCTION 



This article examines the motion of stellar systems in "Einstein-aether" theory — an al- 
ternative theory of gravity that permits breaking of Lorentz symmetry through a dynamical 
preferred frame. The general theory contains four dimensionless couplings whose values can 
be constrained by comparing the predictions of the theory with observations — in particular, 
observations of binary pulsar systems. It will be demonstrated that all current tests will 
be passed by a one-parameter family whose couplings are "small-enough" — that is, on the 
order of 0.01 or less. Identifying whether there is a viable extension of this family to large 
coupling values requires additional work beyond this article. 

The question of whether the physical world is exactly Lorentz invariant has received 
increasing attention in recent years. This interest is sourced largely by hints of Lorentz vio- 
lation ■„ popular candidates for theor.es of quantum gravty-for .nstance. string theory Q, 
loop quantum gravity [2], and noncommutative field theory [3j. More broadly, challenging 
the rule of Lorentz symmetry means challenging the fundamentals of all of modern physics, 
and doing that is iust plain exciting. 

The rev.w ^ di.cu.e= a w.de var.ety of theoretical models that feature Lorent. 
symmetry violating effects, and observational searches for violations. So far no conclusive 
sign of Lorentz variance has been identified, and very strong bounds exist on the size of 
couplings for Lorentz- violating effects in standard model extensions 0, The effects of 
Lorentz violation in a gravitational context, however, are not covered by these bounds. 

Einstein-aether theory — or "ae-theory" for short — is a classical metric theory of gravity 
that contains an additional dynamical vector field. The vector field "aether" is constrained 
to be timelike everywhere and of fixed norm. The aether can be thought of as a remnant 
of unknown, Planck-scale, Lorentz-violating physics. It defines a preferred frame, while its 
status as a dynamical field preserves diffeomorphism invariance. The fixed norm, which can 
always be scaled to unity, ensures that the aetherpicks out just a spacetime direction and 
removes instabilities in the unconstrained theory Gj. 

Much of past work on ae-theory has focussed on placing observational bounds on the 
values of the four free parameters c„ appearing in the ae-theory action, Eqn. ([7]). Constraints 
have been derived from the rate of primordial nucleosynthesis [7|, the rate of Cerenkov 
radiation [G], the requirements of stability and energy positivity of linearized wave modes Q], 
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There, an expression for 



and the parameterized post-Newtonian (PN) form of the theory , Q . A summary of 
these constraints was presented in [8!], where it was shown that they are met by a large 
two-parameter subset of the original four-parameter class of theories. 

Additional constraints on the c„ come from observations of binary pulsar systems. Study 
of the predictions of ae-theory for binary pulsars was begun in 
the rate of radiation damping in A^-body systems was derived to lowest non-trivial PN order 
and neglecting effects due to strong fields in the vicinity of the bodies. It was shown that a one- 
parameter subset of the two-parameter family allowed by the collected constraints discussed 
in would pass tests from binary pulsar systems if there were justification for ignoring 
strong field effects. That neglect is dangerous, though, since the fields inside neutron star 
pulsars should be very strong. Justification requires an unclear assumption on the values of 
the On- 

In this article, I will incorporate strong field effects on binary pulsar systems, calculating 
the PN equations of motion and the rate of radiation damping of a system of strongly 
self-gravitating bodies. The effects will be handled via an effective approach in which the 
compact bodies are treated as point particles whose action contains nonstandard couplings 
that depend on the velocity of the particles in the preferred frame. The effective approach 
to A^-body dynamics in relativistic gravity theories has previously been ^employed in pure 
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131, 



general relativity (GR) and other alternative theories; see for example 
The new interactions are parametrized by dimensionless coefficients, or "sensitivities" , whose 
values can be calculated for a given stellar source by matching the effective theory onto the 
exact, perfect fluid theory. Prior work ll| reveals just the form of the "first" sensitivity at 
lowest order in the self-potential of a body. 

The expressions obtained can be used to constrain the allowed class of ae-theories. Obser- 
vations of binary pulsar systems allow for measurement of "post-Keplerian" (PK) parameters 
that describe perturbations of the binary's Keplerian orbit due to relativistic effects. These 
parameters are mostly "quasi-static" ones, whose expressions can be derived from the non- 
radiative parts of the PN forms of the gravitational fields and the effective equations of 
motion for the bodies. In addition, there is the radiation damping rate, whose expression 
depends on the radiative parts of the fields. The ae-theory expressions for the PK param- 
eters differ from those of pure GR in that they depend on the c„, the sensitivities, and 
the center-of-mass velocity of the system of bodies. Stating precise constraints for general 
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c„ values will require work beyond the scope of this article — specifically, what is needed is 
a method for dealing with dependence on the unmeasurable center-of-mass velocity and a 
calculation of the values of the sensitivities of a given source. 

For the time being, a few comments can be made, which will be defended below. A crucial 
piece of information learned by comparing the weak field limit of the effective theory with 
the weak field limit of the perfect fiuid theory ll| is that the sensitivities will be "small". 
That is, they will be at least as small as {G Nfn / d)"^ , where m is the body's mass and d its 
size, times a c„ dependent coefficient that must scale at least linearly with c„ in the small 
Cn limit. For neutron stars in pure GR, (GNfn/d) ~ (0.1 ~ 0.3); it is reasonable to expect 
something similar in ae-theory based on studies of stellar solutions [l8i] and the fact that 
ae-theory is generally "close" in the small-c„ limit, to GR plus a non-dynamical vector field. 

It then follows that bounds on the magnitude of violations of the strong equivalence 
principle 19| constrain the c„ dependent factor to be less than {0.01){G]\fm/d)~'^ . It further 
follows that the strong field corrections fall below the level of current observational uncer- 
tainties when \cn\ ^ 0.01 and the two conditions that match the ae-theory PPN parameters 
to those of GR are imposed. Thus, weak field analysis 11| suffices for small enough c^, and 
implies the existence of a one-parameter family of theories that passes all current tests from 
binary pulsar systems. 

I will now present the strong field formulas. First, the effective particle action is con- 
structed, and the exact field equations are defined in Sec. [Tll The PN expansions of the 
metric and aether fields are then given in Sec. IIII Al and used to express the PN equations 
of motion for a binary system in Sec. IIII B[ The rate of radiation damping is then determined 
in Sec. IIVI Comments on dealing with center-of-mass velocity and sensitivity dependence 
are given in Sec. |Vl along with the argument for the viability of the weakly coupled family 
of theories. 



I follow the conventions of Wald [20|. In particular, I use units in which the fiat space 
speed-of-light c = 1, and I use metric signature (— , +, +, +). This signature is opposite to 
that employed in HI], but it is much more convenient for calculations involving a time- 
space decomposition. The ae-theory action is defined here in such a way as to permit 
easy comparison between ll| and this article. The following shorthand conventions for 
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combinations of the Cn will be used: 



Cl23 = Ci + C2 + C3, 



Cl4 = Ci + C4, 



C± = Ci ± C3. 



(1) 

(2) 
(3) 



When covariant equations are expanded in Minkowskian coordinates, the following con- 
ventions are observed. Spatial indices will be indicated by lowercase Latin letters from the 
middle of the alphabet: k, . . . . One exception is when the coefficients Ci^2,3,4 are referred 
to collectively as c„. Indices will be raised and lowered with the flat metric rjab- Repeated 
spatial indices will be summed over, regardless of vertical position: Ta — X^j^^ ^Ta. Time 
indices will be indicated by a 0; time derivatives will be denoted by an overdot: / = dof- 

II. EFFECTIVE ACTION AND FIELD EQUATIONS 
A. Particle action 

The aim of this work is to treat within ae-theory a system of compact bodies that poten- 
tially possess strong internal gravitational fields. The complicated internal workings of the 
bodies will be dealt with via an effective approach that reproduces the bulk motion of the 
bodies and the fields far from them. Each body will be treated as a point particle with the 
composition dependent effects encapsulated in nonstandard couplings in the particle action. 

The form of the effective action can be deduced from the following considerations. The 
one-particle action Sa will have the rough form = —'ni J dtO, where the integral is along 
the particle worldline parametrized by t, rh has dimensions of mass, and (9 is a sum of 
dimensionless local scalar quantities. The fundamental theory has only one dimensionful 
parameter G. For a first approximation, the spin of the body can be neglected. Derivative 
couplings in the particle theory are then suppressed by powers of (d/R), where d is the 
size of the underlying finite-sized body and R is the radius of curvature of the background 
spacetime. In addition, Sa presumably reduces to the standard free particle action if the 
particle is comoving with the local aether and must be invariant under reparametrization of 
the particle worldline. 
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These considerations imply the following one-particle action: 

SA = -mA j dTA{l+aA{u%a + l) + '^{u%a + lf + ■■■). (4) 

where A labels the body, ta is the proper time along the body's curve, f " is the body's unit 
four-velocity, and u"" is the aether. The quantity u°'Va expressed in a PN expansion with the 
aether purely timelike at lowest order, is of order f^, the square of the velocity of the body 
in the aether frame. By assumption, v"^ is first PN order (IPN). The IPN corrections to 
Newtonian equations of motion will follow from the part of the action that is tua x (2PN), 
so only the terms in Sa written above are needed for current purposes. For a system of 
particles, the action is given by the sum of copies of Sa- 

This action can be thought of as a Taylor expansion of the standard worldline action, 
but with a mass that is a function of 7 = —u°'Va'- 

Sa = - dTrfiAll]- (5) 



The expansion is made about 7 = 1. The parameters a, a' are then defined as 

aA = -^j^|.=i, aA = aA + aA + aA. = d(h^l-=^- 

This form of Sa suggests that that a a, o'a can be determined by considering asymptotic 
properties of perturbations of static stellar solutions. 



B. Field equations 

The full action is the four-parameter ae-theory action S 

S = j d'x (i? - K-'.^Vau'^Vbu'' + A(m"m V + 1)) , (7) 

plus the sum of copies of S*^ @, retaining only the terms explicitly written above. Here, 

K"',, = {cig'^'gcd + C2S:6', + C3(5^5,' - c.u'^u'g.a) • (8) 
While the sign of the C4 term looks awkward, it permits easier comparison with the results 

of [nl. 



The field equations are then as follows. There are the Einstein equations 

Gab - Sab = ^T^GTab, (9) 



where 



Gab — Rab — -^Rgab, (lO) 

Sab =Vc{K^^''ub) - K^^a^h-j - K^ab)U'') 

+ Ci ( VaMcVfeW'^ - VcU^V^fc) + Ci{u''\/cUa){u'^VdUb) (ll) 
+ XUaUb + ^gab{K''a^cU'^), 



with 



K\ = K'^\,VbU^ (12) 

and T"^ is the particle stress tensor 

T'^' = ^-4^A [A\v-av\ + 2A\u^-vl] , (13) 

A 

with a covariant delta-function 

= (14) 

and 

A\ = l + aA-^{Mr-l), (15) 
^i = -o^A-o^;i(M,t;^ + l). (16) 

The aether field equation is 

VbK'"' = Ci{u'VcUb)V^u^ + Xu" + 87^G'cT^ (17) 

where 



= J]m^5^A>^. (18) 

A 

Varying A gives the constraint QabU^v!' = -1. Eqn. ^ can be used to eliminate A, giving 

A = -u''{VbK\ - c^{Vav'){u'V,Ub) - STrGa,). (19) 

The covariant equation of motion for a single particle has the form 

VbTf - Vb{{(yATu') - ((T^)fcV«n'' = 0, (20) 

where and {ctaY are the one-particle summands in (fT3l) and ( fTSl) . This can be written 
more explicitly as 

v\Vb{A\v\ + A>'^) - A>AbV"u'' = 0. (21) 



III. POST-NEWTONIAN EXPANSION 



A. Fields 



The PN expansion of the fields can be determined by iteratively solving the field equations 
in a weak field, slow motion approximation [s], [2l|. A background of a fiat metric and 
constant aether is assumed, and a Lorentzian coordinate system with the time direction 
defined by the background aether is chosen. Following the procedures of [8] gives 

GNrriA „ G%mAmB ^ \- G%mArhB 



= - 1 + 2 y - 2 y ^N"^^"'^ - 2 y 

+ 3 2^ — : — ^^(1 + 



^ A,B ^ " A,B^A ^ ^" 



ta 

9ij =(l + 2X1^77^)'^*J' 



(22) 



E GnTTLa i , T->+GNmAf j j i 

Ba— < + 2^ B\—^{vyj,)r\, 

A A ^A 



where r\ = — x\, r^^ = x\ — x^. 



Bi = ±1 ± - 2a,) (l + ^|^..) " > + a,)(l + i^.. ). (23) 



and 



= G, (24) 

Z — Ci4 

2ci — c+c_ 

^ _ oi _ (ci + 2C3 - C4)(2ci + 3C2 + C3 + C4) ^2g-) 
2 (2 — Ci4)Ci23 

The numerical values of the PPN parameters ai and 02 are constrained to be very small 
by weak field experiments, yia analysis that allows for a possible lack of Lorentz symmetry 



in the underlying theory 



2l| . There are two independent pairs of conditions on the c„ that 



will set tti and 02 to zero. One pair is 



c. = _M±£l£i:i£l, = (27) 

3ci ci 

The other is c+ = C14 = 0. With this second pair, the spin-1 and spin-0 wave speeds diverge 
(Sec. [IV|l : also, the spin-0 linearized energy density vanishes while that of spin-1 remains 
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finite 



22|. Observational signatures of tliis beliavior have not been worked out, and I will 



not consider these conditions further here. Hence, the first pair of conditions is assumed 
below whenever attention is restricted to the case of vanishing ai and 02 • 
The aether to order of interest is 



" =1+ > , 

— r A 



(28) 



where 

= '-^^'^ ^ -^l (c, + 2c3-c,) + ^2^) 

The results of this section are equivalent to the weak field expressions obtained in fs] when 
a A is set to zero. 



B. Post-Newtonian equations of motion 

The equations of motion for the system of compact bodies follow by expressing the exact 
result (120|) in a PN expansion using the forms of the fields given above. The Newtonian order 
result can be used to define the effective two-body coupling Q and the "active" gravitational 
mass m: 

■ i _ —GN'fnB i _ —Qab'itT'B i 

U + (^A)rAB B^A ^AB 

with the two-body coupling 

Sab = 7:— TT^- r, (31) 

(1 + aA){l + (Tb) 

and the active gravitational mass 

niB = {1 + crB)mB- (32) 

These definitions arise by requiring that Qab = Gba and that mB/rfiB depend on just (Tb- 

Using the Newtonian result and continuing with the expansion leads to the IPN equations 
of motion, expressed here just for the case of a binary system: 
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Gm2 



r \ 1 + o"2 / r 



2V ' 1 + 
-2Dv{vi + ?,{E - D){v{pf 



(33) 



i (vir^ (4 + 3ai - - 3(1 + ai)vi 



+ ?;^(2L)t;jf^' - 2Et;^F) 



where Q = Q12, r"^ = r\ — and 

Z? = -i(8 + ai)(l 



c_ , , (1 — c_) 

— (CTI + 02) H criO-2 



-E = -- - -(ai - 2a2) fl + / 

2 4^ '^V ^ (C1 + 2C3-C4) 



(0-1 + 0-2) - 



2c 



0-10-2 



123 



(34) 
(35) 



The expression for v\ is obtained by exchanging all body-1 quantities and body-2 quantities, 
including the switch r* — > — r*. 



The "Einstein-Infeld-Hoffman" Lagrangian 12| — that is, the effective Lagrangian ex- 
pressed purely in terms of particle quantities — can be determined by working backwards 
from the equations of motion. It is 

L = - (mi + 777,2) + ^("^1^1 + m2vl) 

^ ( ( . o\ \ , /, (To 



+ 1 



1 + Ori 



1 - 



1 + (T2 



+ 



Qm\m2 
r 

1 ( Qmx 
2V r 



l + ^(^(l + (Ti)t;2 + (l + a2)w 
l + .2) + ^(l + ai; 



(36) 



This La gra ngian is not Lorentz invariant unless aA = o'a — 0- This follows from the analysis 



of Will 



231 ] and the list of criteria therein. In particular, the action and the equations of 



motion depend on the velocity of the system's center of mass in the aether frame. 



IV. RADIATION DAMPING RATE 



The radiation damping rate is the rate at which the particle system loses energy via 
gravity-aether radiation. This energy loss manifests as a change in the orbital period of a 
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binary system, equating the energy radiated to minus the change in mechanical energy. The 
expression for the rate in the effective particle theory can be determined by adapting the 
methods of which were used to find the rate for a system of weakly self-gravitating 



perfect fluid bodies in ae-theory. It will be convenient to introduce the parameter sa 



sa = 0-^/(1 + o-a), 



(37) 



and to work with the active gravitational mass 



ruA 



[1 + (TA)rhA = ruA/ (1 - Sa)- 



(38) 



A. Wave forms 



The method of begins by assuming a background of a fiat metric and constant aether, 
with a coordinate system with respect to which the background metric is the Minkowski 
metric T]ab and the background aether is aligned with the time direction. The metric and 
aether perturbations are then decomposed into irreducible transverse and longitudinal pieces. 
The spatial vectors and /loi are written as: 

hoi = 7i + 1,1 = + (39) 

with 7j^j = i'\ = 0. The spatial metric hij is decomposed into a transverse, trace-free tensor, 
a transverse vector, and two scalar quantities giving the transverse and longitudinal traces: 

hij = (pij + ^Pijlf] + 20(ij) + (40) 



where 



(41) 



and 

P^Af]=^^3f,kk- f,f, (42) 

hence, Pij[f]j = 0, and ha = (/ + 0),^. Further, define 

F = In- (43) 

The list of variables then consists of a transverse-traceless spin-2 tensor (pij, transverse spin- 
1 vectors 7j,i/*,0j, and spin-0 scalars 7, i/, F, 0, /ioo, and The Lorentz gauge, or any 
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obvious extension of it, does not usefully simplify the ae-theory field equations. Instead, the 
following convenient conditions will be imposed: 



= M*. = hoi^i = hi[j^k]i, (44) 

or equivalently, 

= z/ = 7 = 0i- (45) 
Because (pi is transverse, these constitute just four conditions. 



Following the field equations can then be linearized and expressed in terms of the 
above variables, and sorted to obtain a set of wave equations with matter terms and nonlinear 
terms as sources. Having done this, the linear contributions can be seen by inspection to 
satisfy a conservation law. This fact implies the existence of a conserved source r'^'' 

^ab ^ j.ab _ ^a^b ^ ~ab ^ ^^g^ 

where T"^ and a" are as defined in Eqns. f|T3l) and f|T8|) . and f'^'' is constructed from nonlinear 
terms — its precise form will not be needed. The non-symmetric r"* satisfies the conservation 
law with respect to the right-index only: r"^^ = 0. The corresponding conserved total energy 
E and momentum P* to lowest PN order are 

E= fd'xT'' = J2^A = J2il- SA)mA, (47) 

A A 

= f d^xT'^ = ^mAv\. (48) 

Conservation of P* means that the system center-of-mass X* defined via rriA 

= (49) 

22a^a 

is unaccelerated to lowest order. 

The field equations reduce to the following. For spin-2, 

1 ■ 



w. 



2 y^ij (Pij,kk 



IGnGr^'^, (50) 



2 



where TT signifies the transverse, trace-free components, and 

1 



wl = (51) 
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For spin-1, 



+ 7i) = ./^""f^ (c+^o + (1 - c+)aO^, (52) 
(c+z^^ + 7^),fcfc = -167rG'rT, (53) 



where T signifies the transverse components, and 



2 2ci — c+c_ , . 

= — r — . (54) 

' 2(1-C+)C14 ^ ^ 



For the spin-0 variables, the constraint gives to hnear order 



M° = 1 + hoo. (55) 



Non-hnear corrections to this are of uninteresting order, as explained in more detail in 11 1. 
The other equations are 

—F - Fkk = [rkk r^fc + — Too) , (56) 

Wq 2 - Ci4 Ci23 Ci4 

{F - Ci4/ioo),fcfc = -167rGroo, (57) 
(1 + C2)Fi + ci230,fcfci = -167rGr^o, (58) 

where L signifies the longitudinal component, and 

(2 — Ci4)Ci23 



(59) 



(2 + 3C2 + c+)(l - C+)ci4 

All these equations can be solved formally via Greens function methods, and the resulting 
integrals expanded in a far field, slow motion approximation. The expressions can be further 
simplified using the conservation of r"^. A result that holds within the approximation scheme 
is that for a field ip satisfying a wave equation with speed w evaluated at field point = |x|n* 
with only outgoing waves, 

wip^i = -Tpfi'. (60) 

Also, differentially transverse becomes equivalent to geometrically transverse to n\ 

The results to lowest PN order and ignoring static contributions are as follows. For spin-2, 

= (61) 
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where the right-hand side is evaluated at time (t — \x\lw2) and the quadrupole moment Qij 
is the trace-free part of the system's second mass moment lij: 



(62) 



For spin-1 variables, 



-2G 



1 



V 



\x\ 2ci — C+C_ Vwi 1 — C+ 

7i = -c+v\ 



Qi3 + - 



(63) 
(64) 



where the right-hand side of the first equation is evaluated at time (t — |x|/wi), Qij is the 
trace- free part of the rescaled mass moment Xif 



(65) 



and 



For spin-0 variables, 



-2G Ci4 

\x\ 2 — Ci4 



A 



2^2 — ttl „\ ^ 2 /S 



2{2c+ — C14) 3woCi4 



^00 — — 

Ci4 
1 +C2 



C123 



(66) 



(67) 

(68) 
(69) 



where the right-hand side of the first equation is evaluated at time {t — \x\/wq), and / = la, 

At this point, the expected smallness of the sensitivities, mentioned in the introduction, 
can be explained. One should take the weak field limit {sa —>■ "small") of the above wave 
forms and compare them with the perfect-fluid theory wave forms determined in ll|. The 
only Syi-dependence that remains at potentially leading order is in Sj. Comparing fl66|) with 
Eqn. (85) of ll| indicates that in the small sa limit. 



= ("1 - ^"2) hO{- 

6 TTLa 



d 



(70) 



14 



where Qa is the binding energy of the body: Q/m {GNfn/d), where d is the characteristic 
size of the body. The imphcation is that when ai = 0:2 = 0, s must scale as (GArm/d)^, 
times a c„ dependent coefficient. This coefficient should scale at least linearly in c„, in the 
Cn limit, to ensure finiteness of the perturbations. 



B. Damping rate expression 



For the next step, the wave forms are inserted into an expression for the rate of change 
of energy S. This expression can be derived via the Noether charge method of Iyer and 



Wald 



24 



251], using the ae-theory Noether charges derived in 26|, with the result 
-1 /■ 1 



^ (2ci-c+c_)(l-c+) ^ .^ . ^ 2 - ci4 



Ff) + O, (71) 



16nG J V 2tU2 '""^ '"^ ' Wi ' ' ■ 4woCi4 J 

where O is a total time- derivative that will be argued away in a moment. 

Using the above results for the wave forms, performing the angular integral, and ignoring 
O gives 
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^ = -G^(^Q,,Q,, + -fQ^,Q^J + ^Q^,Q^, + 1 1 + 1 1 + B^I I + CE^E'), (72) 



where 



A 



Ci4 



2 \W2 (2ci — c^c^y Wi 6(2 — C14) V 2(20+ — C14)/ Wq 



2a2 — tti \ 2 1 



A, 



2^2 — \ 1 



(2-Ci4)c+ 1 ^ / ^ ^ 

2ci — c+c_ Iff V Q{2c^ — C14) J Wq 

1 /2-C14 1 11 



3 / ' 



Ci4 



^ _ Ci4 / 2^2 - y 1 

^ ~ 72 V2(2c+ - Cm)/ wq'' 



+ — + A S+ sr^ "M —1, (73) 

(74) 

A, = — ( ^— - ), (75) 

(76) 

^ 2a2 - ai 1 
' 12(2c+ - C14) w^o" 

= (78) 

6Ci4 

C = -^(iz£M+M. (79) 

3ci4 V wg/ 

The coefficients ^1, Bi, and C are respectively identical to A, B, and C of llj. Taking 
the weak field limit corresponds to retaining only the Ai, Bi, and C terms, and invoking 
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the relation (ITOll for sa in the C term. In the case that ai = 02 = 0, Bi vanishes, as does 
sa in the weak field limit. The weak field damping rate in this case then contains only a 
quadrupole contribution and is identical to the GR rate when Ai = 1. This remaining curve 
of Cn values intersects the range of values allowed by collected constraints considered in j^, 
as illustrated in Figure 1. Thus, this curve gives a one-parameter family of viable ae-theories 
if the weak field results alone are sufficient. 




0.2 0.4 0.6 0.8 



c + 

FIG. 1: Class of allowed ae-theories, if strong field effects in binary pulsar systems can be ignored. 
The four-dimensional c„ space has been restricted to the (c+,c_) plane by setting the PPN pa- 
rameters ai and 02 to zero via the conditions (j27p . The shaded region is the region allowed by 
primordial nucleosynthesis, Cerenkov radiation, linearized stability and energy positivity, and PPN 

n 

constraints, demarcated in |^]. The dashed cm've is the curve along which binary pulsar tests will 
be satisfied, assuming ae-theory weak field expressions. Specifically, it is the curve along which 
^1 = 1 in the ai = 02 = case, so that the damping rate (j72]) is identical to the quadrupole for- 
mula of general relativity. Along both this curve and the boundary of the allowed region, c_ — > 00 
as c+ — > 1. The curve remains within the allowed region for all c+ between and 1. As explained 
in Sec. |Vl strong field effects may lead to system dependent corrections to the binary pulsar curve 
for large c^; however, all such curves will coincide with the weak field curve for |c„| < 0.01 given 
current observational uncertainties. 
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To simplify the expression (172|) , it is crucial to note that the damping rate is calculated to 
lowest PN order using the Newtonian results for the motion of the system. Thus, the system's 
motion can be decomposed into a uniform center-of-mass motion — recall the conservation 
of P* — and a fixed Keplerian orbit in the center-of-mass frame. Since the motion is steady- 
state, the damping rate must have no secular time dependence. This observation implies 
that secular terms in £ arising from Xij, see below, must cancel with secular terms in O. 
In addition, the non-secular portion of O will average to zero when a time average of the 
damping rate over an orbital period is taken, since it is a total time derivative. Thus, O can 
be discarded. 

Hence restricting attention to a binary system, and taking a time average over an orbital 
period, the expression reduces as follows. First, define the quantities 

m = mi + m2, fiA = mA/m, fi = mim2/m, (80) 

and the vectors 

= x{ - 4, v' = r\ (81) 

X' = ^ix\ + /i24, V' =X\ (82) 

To Newtonian order, = —{Qm/r'^)f'^, and = 0. /y can be diagonlized 

lij = l^r'r^ + mX'X\ (83) 

so that 

J.. = '^^JfUl^Sr'f^r - Av^^')). (84) 

As for lij, 

= fi{sifi2 + S2fii)rV + m{sifi, + S2fi2)X'X^ + 2fi{si - S2)r(*X^), (85) 

and 

'iij = S'i'ij - GV^'t^^ + 2^(si - S2yr^'X^\ (86) 

where 

5 = Si/i2 + S2/il, (87) 

and 

E. = (si-.2)^r\ (88) 
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Terms in Xjj with X* dependence are secular; following the discussion above, they can be 
discarded. 

Substituting into Eqn. (1721) and imposing the time average gives the final expression 



£ = -G 



N 



X 



— (^1 + SA2 + S^A3){12v^ - llr') 
15 

+ 4(Si + SB2 + s^B^y 

6 



+ {si - S2) (^(A + 2SA3){3v'V' - 2V¥v^f') + 12{B2 + 2SB3)V'fVf'^j ), (89) 
where the angular brackets denote the time average. 



V. OBSERVATIONAL CONSTRAINTS 



A. Center- of- mass velocity dependence 



While the aether frame center-of-mass velocity of a binary system is not directly 
measurable, dependence of a binary systems's motion on should actually be beneficial 
for constraining the theory. This is because constraints arise from a failure to observe 
V" dependent effects. It may be possible to formulate such constraints without having to 
determine the physical frame, as in the manner of bounds on the PPN parameter 02 • The 
presence of alignment between the sun's spin axis and the ecliptic plane signals the absence 
of frame dependent effects, and leads to a strong bound of |a2| < 4 x 10^'' 2l|. This 
argument does require the assumption that the component of the preferred frame in the 
sun's rest frame is not conveniently aligned with the sun's spin axis; such an assumption 
may generally be required for similar arguments. For example, dependence should cause 
a binary's orbital plane to precess, but not if happens to be normal to the plane. 

An assumption on the order of magnitude of the norm V is necessary to justify the use 
of just the leading PN order expressions for the PK parameters when applied to observed 
binary systems. The validity of the IPN expressions depends on whether corrections of 
relative order and {V^/v"^) are smaller than observational uncertainties. Terms of order 



are negligible for all observed systems, for now, although the "double pulsar" [27| is 
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pushing this hmit. For all but the double pulsar, v"^ ~ 10~^, and uncertainties are at least 
a thousand times this [19]. The double pulsar PSR J0737-3039A/B is the so-far unique 
binary containing two pulsars. The orbital velocity is high, v"^ ~ 10~^, and the presence of 
two pulsars happens to make measurement of system parameters much easier and thus more 
precise — the smallest relative uncertainty is 10~^ on the rate of periastron advance. The v"^ 
corrections are therefore small enough for now, but it is expected that precision will increase 



19|. 



to probe the next PN order within the next 10-20 years 

The dependent terms must feature c„ dependent factors, since it is known that there 
is no center-of-mass velocity dependence at next PN order in pure GR 2l|]. Ignoring those 
factors for the moment, validity of leading PN order for the double pulsar requires that 
{V^/v"^) < 10"^, giving V"^ < 10-^■^ or (V^/v^) < 10°'^ ^ 3. For other systems, given 
uncertainties ranging from (10~^ ~ 10~^), the conditions are {V^/v^) < (10~^ ~ 10~^), 
giving < (10-3.5 _ 10-2-5)^ or (V'^/v'^) < (lO^-^ ~ lO^-^) ^ (300 ~ 30). Presumably, the 
Cn dependent factor actually goes to zero as some positive power of c„, so V can be larger 
in the small c„ limit. A reasonable first guess for the aether frame is the rest frame of the 
cosmic microwave background. A typical velocity for compact objects in our galaxy in this 
frame is V'^ ~ 10~^, so the restriction on V is met. 



B. Constraints in the small coupling regime 

A formula for the sensitivities for a given source should be obtainable by comparing the 
strong field results of this article with analogous results in the exact perfect fluid theory. 
Higher order terms in the exact theory must be calculated, though, since the leading order 



results of 



ll| only give the 0{Gnt^/ d) part of s expressed in f lTOl) . The calculation can 



be done in the case of a single body that is static except for a constant aether frame 
velocity, by, for example, continuing the iterative procedure used to determine the PPN 
parameters js, 23]. The process may be lengthy, but straightforward. 

I have shown that the sensitivity of a body will scale with the body's self-potential like 
[3[cj^{G jqT^ / d)"^ , where (3 is some c„-dependent coefficient that scales at least as fast as c„ in 
the small c„ limit. Even in the absence of a formula for the sensitivities and precise knowledge 
of center-of-mass velocities, two useful comments can be derived. First, a constraint can be 
roughly stated: \(3\ < (0.1 ~ 1). Second, there exists a one-parameter family of theories that 
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passes all current constraints, obtained by restricting to c„ with magnitude less than roughly 
0.1 and imposing the two PPN conditions, the one weak field damping rate condition, and 
the collected non-binary conditions. 



The condition that \j3\ < (0.1 ~ 1) follows from constraints 19| on the magnitude of vio- 
lations of the strong equivalence principle — that is, that a body's acceleration is independent 
of its composition. A violation would lead to a polarization of the orbit of pulsar systems 
due to unequal acceleration of the binary bodies in the gravitational field of the galaxy. The 
observed lack of polarization in neutron star-white dwarf systems leads to a constraint that 
can be stated here as s < 0.01, where here s is the sensitivity of the neutron star in the 
considered pulsars. Assuming that {0^171/ d) (0.1 ~ 0.3) for the pulsar, as it is in GR, 
the constraint on the size of (3 arises. It is possible that when the weak field conditions are 
imposed, (3 will automatically satisfy the above inequality; certainly it will in the small c„ 
regime when |c„| < 0.01. 

The statement, that current tests will be satisfied if the weak field conditions are imposed 
and the remaining degree of c„ freedom satisfies |c„| < 0.01, can be derived by considering 
the battery of binary pulsar tests. First, consider tests that probe only the quasi-static 
PK parameters — that is, all but the damping rate. The tightest quasi-static test comes 
from the double pulsar [27]. The relative size of the strong field corrections to the weak field 
expressions will be O(s^), while the prediction of GR has been confirmed to within a relative 
observational uncertainty of 0.05%. Requiring s < 10^^ and assuming that (GNm/d) ^ 
(.1 ~ .3) for the pulsars, the condition \cn\ ^ 0.01 arises. Given this and the two conditions 
that set the PPN parameters ai and 02 to zero, all current quasi-static tests will be passed. 

Tests that incorporate the damping rate will also be satisfied by the small-c„ condition and 



the weak field conditions. I note first that for systems in which the damping rate is 



Drobed, 



19 



2ll- 



uncertainty on its measurement dominates uncertainties on quasi-static parameters 
Thus, it is conventional to use the measurements of the quasi-static parameters to solve 
for the mass values of the binary bodies. When ai = 02 = 0, and |c„| < 0.01, so that 
the expressions for the quasi-static parameters are close to those of GR, the predicted mass 
values will also be close. 

Now, the dipole contribution to S can be significant in asymmetric systems where the 
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sensitivity of one body is much larger than the other. The dipole contribution is 

Empoie = -Gj,{{^f)C{s, - S2)\ (90) 

which is of order (Cs^/lOf ^) compared to the quadrupole and monopole contribution, where 
s is the dominant sensitivity. An applicable system is a neutron star-white dwarf binary, 
since for a typical white dwarf, (GNiTi/d) ~ 10"'^. Constraints have been derived 19|] on the 
magnitude of dipole radiation from neutron star-white dwarf binaries PSR B0655+64 and 
PSR J1012+5307 by requiring that the dipole radiation rate be no larger than the observed 
rate. The analysis applied here leads to the condition Cs"^ < 10~^, where s is the sensitivity 
of the neutron star. In the small c„ regime, this translates again to the condition \cn\ ^ 0.01. 

For double neutron star binaries, the dipole rate is further suppressed by the similarities 
of the sensitivities, and the quadrupole and monopole contributions become dominant. The 
tightest test involving radiation is associated with the Hulse-Taylor binary PSR1913+16, 
with a relative uncertainty of 0.2% 19|, |21|. In the small c„ regime, the condition = 1 
matches the leading order damping rate to that of GR. The strong field corrections are of 
relative order s; to be smaller than the uncertainty again requires |c„| < 0.01. 

This upper limit on |c„| will decrease as observational uncertainties decrease. The 
most promising candidate for lowering the limit is the double pulsar: 2PN-order and spin- 
dependent effects should be observable within the next ten or twenty years 271]. Another 
type of system, yet undetected, for which high levels of accuracy could be obtained is a 
neutron star-black hole binary, as the structureless black hole would decrease noise due to 
finite-size effects and mass transfer between the bodies. 

For \cn\ > 0.01, strong field contributions to the expressions for the PK parameters may 
be significant. Those contributions for a given source cannot yet be calculated, so the theory 
cannot be checked against observations. Thus, there is no conclusion yet on the viability of 
large c„ values. If it were possible to calculate precise predictions for a given binary system, 
then each observed system would imply an extension from small to large c„ of the curve of 
allowed values. The only physically viable values would be those for which the curves for 
all observed systems overlapped within error. 
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